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Chiral perturbation theory for lattice QCD at 0(a 2 )* 
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The 0(a 2 ) contributions to the chiral effective Lagrangian for lattice QCD with Wilson fermions are constructed. 
The results are generalized to partially quenched QCD with Wilson fermions as well as to the "mixed" lattice 
theory with Wilson sea quarks and Ginsparg- Wilson valence quarks. 



1. Introduction 

Chiral perturbation theory (%PT) is widely 
used in analyzing lattice QCD data. Since xPT 
describes continuum QCD one has to take the 
continuum limit first before fitting the lattice 
data to the analytic predictions of xPT- Unfor- 
tunately, lattice data for various lattice spacings 
is often not available in order to perform a re- 
liable continuum extrapolation. Nevertheless, it 
is common practice to directly fit lattice data ob- 
tained for non- vanishing lattice spacing to contin- 
uum xPT, thus introducing a systematic error. 

This uncertainty can be analytically controlled 
by including the discretization effects stemming 
from a non-zero lattice spacing in xPT- For ex- 
ample, the 0(a 2 ) taste violations in staggered 
fermion simulations have been systematically ac- 
counted for that way £Q. Fits of MILC simula- 
tion data to xPT including 0(a 2 ) taste violat- 
ing terms show a significantly improved \ 2 0- 
In fact, good fits are not possible without these 
terms. 

Similarly, the leading effects linear in a for sim- 
ulations with Wilson fermions have been included 
in xPT P] and the resulting expressions for the 
pseudo scalar mass and decay constant have al- 
ready been used in fits to simulation data 0]. 
However, since the lattice spacings in current un- 
quenched simulations with Wilson fermions are 
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usually not very small, it is important to extend 
the results in [3] by including the 0(a 2 ) contri- 
butions in xPT 

2. for Wilson fermions at 0(a 2 ) 

The strategy for including the cut-off artifacts 
in xPT is essentially a two-step matching to ef- 
fective theories. We first match the lattice theory 
to Symanzik's effective theory [7| , which is a con- 
tinuum field theory. This effective theory makes 
the cut-off dependence of the underlying lattice 
theory explicit in terms of higher dimensional op- 
erators in the effective action and the effective 
operators. In a second step we construct the chi- 
ral effective theory for Symanzik's effective the- 
ory, including the extra terms due to the non-zero 
lattice spacing. The result is a chiral expansion 
with explicit dependence on the lattice spacing a. 

In order to include 0(a 2 ) effects in xPT we 
need Symanzik's effective action for lattice QCD 
with Wilson fermions through 0(a 2 ). At this or- 
der there are fifteen dimension-6 operators which 
are compatible with the symmetries of the un- 
derlying lattice theory [Sj. The way these opera- 
tors affect the chiral Lagrangian depends on their 
transformation properties under chiral rotations, 
and one can basically distinguish three types of 
operators: 

1. the 0(4) symmetry breaking term 
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2. chiral symmetry conserving terms 
e.g. ^ an d {^l^i^l^) 

3. chiral symmetry breaking terms 
e.g. (ipj 5 i/j)(ipj 5 ip) 

The first operator, which breaks the full 0(4) ro- 
tation symmetry, is not forbidden since the un- 
derlying lattice theory is invariant only under the 
hypercubic subgroup of 0(4). 

The chiral Lagrangian is expanded in powers 
of momenta p 2 , quark masses m and the lat- 
tice spacing a. Generalizing the standard chi- 
ral power counting, the leading order (LO) chi- 
ral Lagrangian contains the terms of 0(p 2 , m, a), 
while the terms of 0(p ,p 2 m,p 2 a 1 m 2 , ma, a 2 ) are 
of next-to-leading order (NLO). Employing this 
power counting the first two types of operators 
given above do not contribute at NLO. 

The first operator, for instance, gives rise to 
terms in the chiral Lagrangian which break the 
0(4) symmetry and are invariant under the hy- 
percubic subgroup only. These terms necessarily 
contain four derivatives and, since they are also 
multiplied by a 2 , are therefore at least of 0(p 4 a 2 ). 
These terms are beyond NLO and can be safely 
neglected. 

The second type of operators, the chiral sym- 
metry conserving terms, modify the coefficients 
in front of already existing terms in the chi- 
ral Lagrangian. For example, the LO coefficient 
f 2 , multiplying the kinetic term tr be- 
comes / 2 (a 2 ) = f 2 + a 2 K + . . ., with some other 
(unknown) coefficient K. This leads to the ex- 
tra term a 2 Ktr c^EdE^ in the chiral Lagrangian, 
which is of 0(a 2 p 2 ) and thus beyond NLO. This 
example can be generalized to all chiral symme- 
try conserving terms and we conclude that these 
terms do not affect the NLO chiral Lagrangian. 

Only some of the chiral symmetry breaking 
terms give rise to new operators in the NLO chi- 
ral Lagrangian. The way they influence the chi- 
ral Lagrangian is obtained by a standard spurion 
analysis and the final result for the contributions 
proportional to a 2 reads 

C[a 2 ] = -a 2 W^(£t + £) 2 -a 2 ^(l]t-E) 2 
-o 2 W^(S t E t + EE). (1) 



Here £ is the standard matrix of Nambu- 
Goldstone fields and the angled brackets denote 
traces over the flavor indices. The parameter 
a = aWo absorbs the leading order low-energy 
constant Wq in order to render the new unknown 
low-energy constants W 6 , W 7 , W 8 dimensionless. 
Note that only three new operators appear at 
0(a 2 ) despite the fact that there are many differ- 
ent operators in the underlying Symanzik action 
at this order. 

We are now in the position to compute NLO ex- 
pressions for physical observables. For example, 
the 0(a 2 ) contribution to the flavor non-singlet 
pseudo-scalar mass is given by 

16 i i 

6M 2 PS [a 2 } = — a 2 (N f W 6 + W s ), (2) 

where Nf denotes the number of flavors. Adding 
this contribution to the terms of 0{a) given in 
[3] we obtain the complete NLO result through 
0(a 2 ) for the pseudo-scalar mass. 

3. Partially quenched QCD 

The previous discussion is readily extended to 
partially quenched lattice QCD with different 
masses for the sea and valence Wilson fcrmions. 
Described by a lattice action with sea, valence 
and ghost quarks, the Symanzik action through 
0(a 2 ) for partially quenched lattice QCD is ob- 
tained as in the unquenched case, based on local- 
ity and the symmetries of the lattice theory. The 
result for the 0(a 2 ) contribution to the partially 
quenched chiral Lagrangian has the same form 
as in QJ, with the angled brackets now denoting 
super-traces and the field £ reflecting the larger 
flavor content of partially quenched xPT. 

The discussion of "mixed fermion theories" is 
slightly more complicated. Mixed theories are a 
generalization of partially quenched lattice the- 
ories. In addition to choosing different quark 
masses, the lattice Dirac operator is different in 
the sea and valence sector. The concrete example 
with the Wilson Dirac operator for the sea quarks 
and a Dirac operator satisfying the Ginsparg- 
Wilson relation for the valence quarks has been 
proposed in 

Such a combination has various advantages. 
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As a consequence of the Ginsparg- Wilson rela- 
tion the valence sector exhibits an exact chiral 
symmetry and therefore does not suffer from the 
shortcomings of Wilson fermions due to their ex- 
plicit chiral symmetry breaking (additive mass 
renormalization, cut-off effects linear in a, etc.). 
In particular, simulations are possible with va- 
lence quark masses much smaller than accessi- 
ble with Wilson fermions ^HI- This should al- 
low numerical simulations deeper in in the chiral 
regime of partially quenched QCD. Finally, mixed 
fcrmion simulations offer a cost-effective compro- 
mise towards full unquenched simulations using 
Ginsparg- Wilson fermions, since mixed simula- 
tions only require the computation of correlators 
in the background of already existing unquenched 
configurations generated with Wilson fermions. 

Deriving the 0(a 2 ) corrections in the chiral La- 
grangian for the mixed fermion theory with Wil- 
son sea and Ginsparg- Wilson valence quarks goes 
along the lines outlined above. However, due to 
the use of Ginsparg- Wilson fermions for the va- 
lence quarks there are no G{a) terms and fewer 
0(a 2 ) terms in the valence sector of Symanzik's 
effective action for the mixed theory. As before 
one can show that the 0(4) symmetry breaking 
terms as well as the chiral symmetry conserving 
operators in Symanzik's action do not contribute 
to the NLO chiral Lagrangian. A spurion analy- 
sis for the chiral symmetry breaking terms gives 
the following 0(a 2 ) contribution for the chiral La- 
grangian: 

£[a 2 ] = -a 2 Ws(P s ^ + SP S ) 2 
-a 2 W$(P s tf - SP 5 ) 2 
-a 2 W^{P s ^Ps^ + SP S EP S ) 
-a 2 W M {T 3 ^T 3 ^) (3) 

The operators Ps and Py denote projection oper- 
ators to the sea and valence sector, respectively, 
and T3 equals Ps — Py (as before, a = aWo). 

The first three terms stem entirely from the 
Wilson fermion sea sector. They involve sea sec- 
tor fields only (note the projector Ps) and they 
are the same three operators as in JU . The addi- 
tional operator mixes the sea and valence sector. 
The reason for the presence of the extra term is 
the reduced flavor symmetry group of the mixed 



theory. The use of different Dirac operators for 
sea and valence quarks does not contain transfor- 
mations between sea and valence fields and the 
additional term (t^Ht^T^) is not forbidden by the 
symmetries of the theory. 

The presence of this extra term in the chiral La- 
grangian does not entail that chiral expressions 
for observables in the mixed theory depend on 
more free parameters. For example, taking ac- 
count of the terms in J3J in the calculation of the 
flavor non-singlet pseudo-scalar mass one finds 
that its 0(a 2 ) correction actually vanishes. This 
is not that surprising. Exact chiral symmetry in 
the massless case for Ginsparg- Wilson fermions 
implies that the pseudo-scalar meson mass is pro- 
portional to the quark mass and vanishes in the 
chiral limit. Hence any lattice spacing correction 
to Mp s is suppressed by at least one factor of the 
valence quark mass mvai- The largest lattice cor- 
rection quadratic in the lattice spacing is there- 
fore of O (my aid 2 )- This example illustrates that 
certain beneficial properties of Ginsparg- Wilson 
fermions are preserved even in the presence of a 
"non- Ginsparg- Wilson" sea sector. 
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